Bell-type inequalities for cold heteronuclear molecules 
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We introduce Bell-type inequalities allowing for non-locality and entanglement tests with two 
cold heteronuclear molecules. The proposed inequalities are based on correlations between each 
molecule spatial orientation, an observable which can be experimentally measured with present 
day technology. Orientation measurements are performed on each subsystem at different times. 
These times play the role of the polarizer angles in Bell tests realized with photons. We discuss 
the experimental implementations of the proposed tests, which could also be adapted to other high 
dimensional quantum angular momenta systems. 
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Entanglement and non-locality are two related pecu- 
liarities of quantum theory. Historical debates [l| and the 
advent of quantum information theory, to which entan- 
glement is an essential ingredient, increased the interest 
paid to those concepts considered to be among the main 
traits of quantum mechanics 0] • The experimental real- 
ization of Bell Q and Bell-type inequality tests Q is of 
extreme importance to confirm the validity of quantum 
theory in a loophole free context [j| . Since entanglement 
and non-locality are related (non-locality implies entan- 
glement but not all entangled states are non-local), such 
tests may equally serve as entanglement witnesses @ . 
The interest of such approach also lies in the fact that for 
arbitrary high dimensional systems one does not dispose 
of necessary and sufficient conditions for entanglement 
characterization. However, detecting quantum proper- 
ties of high dimensional systems is of clear interest in 
atomic, molecular and optical physics. 

A number of experimental proposals are aiming to 
throw some light to these questions by means of Bell- 
type tests performed in different continuous or multi- 
dimensional systems: it was shown recently that corre- 
lated atoms originated from the dissociation of a molec- 
ular Bose-Einstein condensate can be used to test the 
Einstein-Podolski-Rosen paradox as originally formu- 
lated [H, 0] ■ The dichotomization of measurement results 
performed in a continuous system can also lead to maxi- 
mal violation of Bell- type inequalities in phase space Q . 
Finally, discrete multi-dimensional systems also allow 
for non-locality tests 0] , and recent experiments show 
the violation of Bell-type inequalities with two effective 
spin-1 particles [l(| ■ In parallel, recent advances in single 
molecule manipulation and detection open the way to the 
controlled creation of molecular entanglement . More- 
over, cold polar molecules confined in optical or mag- 
netic traps have been recognized as a promissing candi- 
date for quantum information processing [l2| , espec ially 
when their rotational levels are used as qubits [Tj, [lj] . 
Rotational states are relatively long lived, allowing for 



short quantum gate implementation times: one can per- 
form about 10 4 gate operations before decohcrcncc takes 
place. This excellent performance, when compared to 
cold collision based quantum gates [l5| , originates from 
the strength of dipolar interactions. 

In the present paper, we propose realistic non-locality 
and entanglement tests for a system composed of two 
cold and trapped heteronuclear diatomic molecules. As 
in usual Bell tests scenarios, measurements are performed 
independently on each molecule by observers placed far 
apart, so that no comunication between them is possible 
during the realization of the protocol. We show in the fol- 
lowing that inequalities built from local realism assump- 
tions can be violated by a set of entangled states using 
measurements of correlations between the two molecules 
spatial orientation. These inequalities can also be used as 
entanglement witnesses [(| when the measurement condi- 
tions are loosen. 

The inequalities derived here follow the original formu- 
lation of Clauser, Home, Shimony and Holt (CHSH) 
For a pair of spin-1 /2 particles or equivalent two-level 
systems, it can be shown using local theories that 
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where a a is the Pauli matrix in the a direction, a and 
a' refer to the first particle while b and b' refer to the 
second one. In order to derive the bound appearing in 
Eq. |T]) one assumes that the measurements performed 
by each observer are independent and their results have 
a probability distribution which is a product of indepen- 
dent probabilities for each subsystem. Such probabilities 
can also depend on some random local variable. Details 
are discussed in many works, as [l6| for instance. It can 
be shown that a collection of entangled states can violate 
Eq. |T]), and this experimental violation was observed us- 
ing photon pairs entangled in polarization [H[ . In this 
case, the directions a, a', b and b' refer to different ori- 
entations of polarizers placed before the detectors. The 
value of the maximal possible violation, 2v / 2, was first 
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derived by Cirel'son [l7| • 

The Bell inequalities we propose here have a structure 
similar to Eq. {!]). They arc based on molecular orien- 
tation correlation measurements instead of spin-like ob- 
scrvables. We show that they are experimentally imple- 
mentable using time delayed measurements. 

Molecular spatial orientation : We suppose that the 
two linear molecules that compose our system behave 
like rigid rotors described by their angular properties 
at a given time t. An arbitrary two-molecule bipartite 
state \ipo) has been created at time t = 0, after which it 
freely evolves. This state can be entangled using tech 
niques as the ones proposed in [l" 
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No interaction 
between molecules is allowed for t > 0. The Hamilto- 
nian describing each individual molecule's free evolution 
is Hi = Jf/h 2 , where i = 1,2. It is expressed in units 
of the rotational energy. Ji is the angular momentum 
operator and therefore the associated evolution operator 
is given by Ui(t) = e-wflit/^ where time t is written 
in units of the rotational period r. Uiit) is therefore 
time-periodic with period 1. The total molecular state 
is thus given by the wavefunction ip(9i, 62, <j>i, 02j t) = 
(6*i, 62, 4>i, (j>2\ip(t)}. 0i and <j>i denote here the polar and 
azimuthal spherical coordinates of both molecules in the 
laboratory frame. For each subsystem of the molecular 
bipartite set, the orientation at time t is defined as the 
average value (4>o\Oi(t)\ipo) of the orientation operator 
Oi(t) = U^it) <g J7 1 _1 (i) cos(0i) l>i(i) <8> U 2 (t). Note that 
here cos(^) is taken as an operator, in perfect analogy 
to the position operator and related functions. Orienta- 
tion as defined above can be experimentally measured by 
recording each molecule's fragments angular distribution 
following a quasi-instantaneous molecular dissociation in- 
duced by ultrafast Coulomb explosion [3] ■ 

With an arbitrary accuracy, each molecule's state (sub- 
scripts will be omitted) can be considered to reside in a 
finite dimensional Hilbcrt space TL generated by the basis 
set {\j,m); j s$ j m ^ x , \m\ < j}, where \j,m) are the 
eigenstates of j 2 and J z . Since it is a priori possible to 
prepare entangled bipartite molecular states with m = 
using ultra-cold atomic photoassociation or Feshbach res- 
onances followed by a dipolar interaction scheme 12, IH , 
we have considered here, for simplicity, that the angular 
momentum projection m is fixed at for each molecule. 
We have checked that our conclusions remain valid even 
if m ^ 0. The dimension of each molecule's Hilbcrt space 
H is then j max + 1, and the states \j,m = 0) will now 
simply be written as \j). The corresponding wavefunc- 
tions are the spherical harmonics (6\j) = Yjo(6). In the 
finite Hilbert space H, the cos 9 operator is character- 
ized by a discrete, non-degenerate spectrum of eigenval- 
ues A„ (0 ^ n ^ jmax), with corresponding eigenvectors 
|A n ), also called orientation eigenstates. The two maxi- 
mally oriented states |+) and |— ) are the two eigenstates 
corresponding to the extreme eigenvalues ±A max , where 



Amax = Max(A„). In the particular case of j max = 1, the 
maximally oriented states can be simply written in the 
basis of the angular momentum eigenstates |0) and |1) as 
|+) oc |0) + |1) and |— ) cx |0) — |1). In this particular case, 
cos 9 oc a x , and we can of course recover the well- know 
results of a two-level system. 

We now describe how orientation measurements per- 
formed locally on each subsystem can be used for non- 
locality tests through inequalities analogous to Eq. fl}. 

Temporal inequalities : The orientation correlation be- 
tween the two molecules is now defined as the average 
value of the two-particle operator C(t\,t2) = Oi(ii) <8> 
02(^2), measured at times t\ and t 2 - Note that the oper- 
ator cos 9 is useful here for entanglement detection since 
it "mixes" different values of j without affecting their 
projection m. Previous works have considered correla- 
tions between different values of the projection m for a 
given (fixed) value of j in other physical contexts 0, [3] . 

The rotation of the molecular axis is described by the 
free evolution operators Uiit). As coherent superposi- 
tions of a finite numbers of \j), the orientation eigenstates 
are not stationnary states of the free Hamiltonian. Time 
evolution therefore modifies these quantum superposi- 
tions, thus changing the molecular orientation, in exact 
analogy to the usual projection of the photon polarization 
with polarizers. We would like to stress here that other 
CHSH inequalities using free evolution instead ofpolar- 
izers were studied in different contexts: in Ref. dEm, 
they allow the detection of entan glem ent between prod- 
ucts of decaying mesons. In Ref. [21|, they reveal quan- 
tum properties of single particles. 

By combining measurements realized at different 
times, one can define, in analogy to Eq. (fT]), the operator 

Bi(t) = C(t, t) + C(t, 0) + C(0, t) - C(0, 0). (2) 

For simplicity, we have used a single time variable t here 
in comparison with the four angles of Eq. (fT]) . We have 
verified that this limitation does not affect the generality 
of our results, since introducing four measurement times 
ti only increases the number of entangled states detected. 
In the framework of a local theory (LT), the operator 
B\ it) obeys an inequality similar to Eq. (JTJ) : 



(Bi(*)> 



LT 



^2A 



Vt > 0. 



(3) 



Without loss of generality, we have also assumed that 
each particle state resides in the same finite dimensional 
Hilbert space TC. We note that Eq. |3]) is valid for all 
possible values of j m ax, and that it can, in particular, 
be extended to the limit j ma x — y +00, in which case 
the spectrum of cos 9 forms a continuum. An interest- 
ing characteristic of the separability threshold ([3]) is its 
dependence on A max . In order to study the spectrum of 
Bi(t), we can numerically diagonalize this operator. We 
obtain, for each time t, its highest eigenvalue, which gives 
the maximal value of (B\(t)), and therefore the maximal 
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FIG. 1: (Color online) Maximal value of (Bi(t)) in the sub- 
space 7i with j ^ jmax as a function of t in units of the ro- 
tational period r. The black, red and blue curves correspond 
to jmax = 1, 2 and 5 respectively. The associated separability 
thresholds are represented by dashed horizontal lines. The 
continuous case associated with the limit j max — ► +oo is also 
shown as a dotted horizontal line. 

violation of Eq. ([3]) . This quantity depends on the dimen- 
sionality of the system, and we compare the amplitude of 
the violation for different values of jmax in Fig.[T] While 
jmax = 1 corresponds to a two-level system, j max = 2 
and jmax = 5 are two examples of a higher dimensional 
system where local realism can, in principle, be violated. 

A first motivating result is that for a broad range 
of times, Max[(Bi(i))] > 2A max , thus violating Eq. ©. 
It proves that (6i(t)) is useful for non-locality tests. 
In particular, in the two dimensional case j max = E 
Max[($i(t))] reaches its maximal value of 2y/2 A max , in 
a similar fashion as for CHSH inequalities with two level 
systems. However, with increasing dimensionality, the 
maximum relative violation decreases. This effect is 
shown in Fig. [I] the maximal relative violation of Eq.j3]) 
decreases from 41% to 29% when j max varies from 1 to 
5. In addition, the highest possible separability threshold 
given by Eq. ^ is obtained in a true infinite dimensional 
space, since in this case A max — > 1. The cos 9 eigenvalues 
then form a continuum, and Eq. ^ becomes a CHSH 
inequality for continuous measurement values. In a real 
experiment, one does sometimes not control the dimen- 
sion of the subspace where entanglement is created, and 
the maximum threshold 2A max — > 2 should then be con- 
sidered. Fig. [T] shows that already for low values of jmax 
it is possible to violate this general threshold. Indeed, for 
jmax = 5, the maximum value of (B\(t)) is 9/4, and for all 
values jmax 5, one can find entangled states violating 
Eq. ([3]). This result can be very useful in the realization 
of an experiment, and it is a consequence of two facts: 
high orientation (high values of A max ) can be obtained in 
reduced angular momentum subspaces [22[ and entangle- 
ment enhances two particle orientation correlations [l3| . 

We have shown that the operator defined by Eq. @ 
allows not only for the realization of Bell-type tests in 



FIG. 2: (Color online) Maximal value of {Bi(t)) as a function 
of t in units of the rotational period r. The solid, dashed and 
dotted curves correspond to j m ax = 1, 2 and 5 respectively. 
The separability and Cirel'son bounds are shown as horizontal 
lines. 

finite angular momentum subspaces, but also when the 
size of the subspace is not a priori known. However, 
the maximal value of the violation decreases with di- 
mensionality. This may render the proposed test more 
difficult in very high dimensions. We can overcome this 
by a dichotomizing procedure, in which a high dimen- 
sional system is transformed into an effective two level 
one. The dichotomization is defined as follows: the 
states | A n ) for which A n = (A„| cos# |A n ) > are said 
to be positively oriented, while those for which A„ 
are considered as negatively oriented. The orientation 
cigenstatcs are thus separated in two classes, |A+) and 
|A_), depending on the sign of their associated eigen- 
value. We now define, for each molecule, the projec- 
tors on the subspaces of positive or negative orienta- 
tion: n± = 53 A |A±)(A±|. The measured observable 
for the molecule i is then H; = n + — n_ . This measure- 
ment consists in counting the asymmetry of the molecular 
angular distribution. For a given single-molecule state 
If) = Ea+ c a+I^+> + Ea_ c aJA-), (<p\IU\<p) can take 
any value in the interval [+1,-1], The total two-molecule 
observable is defined as n = Hi (g) H2 . We refer, as be- 
fore, to two-molecule correlation measurements realized 
at two different times, using H(ti,t2) = Hi(ti) <g> ^(£2) 
where 11, (U) = U" 1 ^)!!.^^). In analogy to Eq. ©, 
we now define the operator 

B 2 (t) - n(t, t) + n(t, o) + n(o, t) - n(o, o). (4) 

Since H(ti,tj) 2 — 1, one can show that the highest 
value {Bjit)) can reach is given by the Cirel'son bound 
2v / 2 [17B ■ In addition, in the framework of a local theory, 
we have 

KB 2 (i)>LTK2, Vi>0. (5) 

Fig. [2] shows the maximum value of (^(i)) as a function 
of time for j max = 1, 2 and 5. For j max = 1 we obtain 
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trivially the same result as with Bi (with a simple scaling 
factor). However, for j max > 1, measuring B2 has some 
advantages over B\ . Cirel'son and locality bounds do not 
depend on dimensionality anymore, remaining valid even 
in the continuous case. Furthermore, Fig.[T] and[2] show 
that the range of time when B2 violates locality is much 
larger than the one of B\. 

We now illustrate the principles of an ideal experiment 
leading to non-locality (and entanglement) tests using 
Eq. ((4|. Our scenario consists of two molecules, created 
in a bipartite state |^o) whose non-local properties one 
wishes to test. The molecules are spatially separated 
and evolve freely so that observers A and B can inde- 
pendentl y m easure their orientation following the prin- 
ciples of[l8(. Measuring the average value of the di- 
chotomized orientation ([¥]) requires only two detectors for 
each molecule, one placed in each hemisphere, determin- 
ing whether molecules are positively or negatively ori- 
ented. Observers A and B perform such measurements 
at different times ti. To ensure non-locality tests are re- 
alized, we impose that the measurements performed by 
each observer cannot causally affect the other's. For that, 
we suppose that A and B are separated by a distance I 
such as I > ct, where c is the speed of light and r the 
molecule's rotational period. For typical values of r, this 
means that the observers should be separated by a dis- 
tance greater than a few hundred /^m. A possible experi- 
mental setup allowing for such molecular separation and 
manipulation are optical tweezers, currently used as indi- 
vidual atoms traps [i^] . The protocol is thus the follow- 
ing: after collecting orientation measurements for differ- 
ent times ti, observers A and B compare their results and 
construct the quantity Using measurements lying in 
the interval determined by a rotational period r, they can 
experimentally infer the violation of (j4]) . In a simpler ap- 
proach, we can also consider the situation where £ < ct. 
In this case, correlations between measurements realized 
at different times and the proposed inequalities can serve 
as entanglement witnesses. This version of the proto- 
col may find direct applications for molecules trapped in 
optical lattices [HI ■ 

In conclusion, we have introduced two experimentally 
realistic distinct high dimensional Bell-type inequalities 
based on the measurement of a quantity which is usually 
considered as being "almost classical" : the molecular ori- 
entation. Both can be violated not only in the case of a 
restricted angular momentum subspace but also in the 
continuous limit. The proposed inequalities also present 
another original feature: orientation measurements are 
realized at different times for each particle. We have 
shown that dichotomization into two classes, namely pos- 
itive and negative orientation, allows for a simple experi- 



mental measurement procedure that can be implemented 
with present day technologies using molecules in optical 
traps. This second approach also increases the number of 
detected entangled states. Our results open the perspec- 
tive of entanglement detection and non-locality tests for 
high angular momentum systems in atomic and molecu- 
lar physics. 
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